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The  bandwidths  for  tnetastable  states  in  a 
simple  model  of  amorphous  semiconductors  is  evaluated 
exactly.  The  region  of  metastability  in  the  space  of 
the  model  parameters  (phase  diagram)  is  deduced  and  the 
subregion  corresponding  to  an  infinite  range  model  is 
indicated.  The  model  predicts  both  paramagnetic  and 
diamagnetic  metastable  states  and  the  relationship 
between  them  is  discussed. 
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In  this  paper  we  further  examine  the  properties  of  metastable 
electronic  states  predicted  for  the  negative  U  model  (Anderson,  1975)  of 
amorphous  semiconductors  by  Licciardello,  Stein  and  Haldane  (1980) 

(see  also,  Licciardello  1980a).  The  presence  of  such  states  have  been  linked 
to  photoinduced  paramagnetism  (Mott,  Davis  and  Street,  1975;  Licciardello 
and  Stein,  1980)  and  conductivity  changes  (Licciardello,  1980b)  observed  in 
a  wide  variety  of  glassy  semiconductors.  Interest  in  the  properties  of 
metastable  states  stems  from  their  defect-like  behavior,  although  the  model 
itself  contains  no  defects  per  se.  Thus  the  theory  provides  an  alternative 
point  of  view  to  the  generally  accepted  defect  picture  of  glassy  materials 
(Kastner,  Adler,  and  Frit2sche,  1976). 

We  consider  a  spatially  random  distribution  of  localized  states  in 
which  the  Hamiltonian  H^  corresponding  to  the  itl1  site  is  written 

Hi  =  Ci"i  +  ubibi  +  U"i+”i4-  “  +  b*)^  (1) 

where  the  quantity  ru  is  the  number  operator  for  an  electron  at  the  i*"*1 
site.  Similarly  n^  (h^)  counts  electrons  at  the  iC^  site  with  spin  up  1 
(down  F)  .  I  have  also  introduced  the  operator  b+  which  creates  a  phonon 
at  the  site  with  frequency  m  .  The  parameters  and  g  are, 

respectively,  the  electron  energy  at  the  i^  site,  the  Hubbard  repulsion 
energy  and  the  electron  phonon  coupling.  We  have  ignored  the  zero  point 
lattice  vibrations  and  employ,  for  simplicity,  an  Einstein  model  for  the 
phonons.  We  also  neglect  the  variation  of  the  coupling  g  and  the  Hubbard 
repulsion  U  from  site  to  site. 

The  Hamiltonian  has  the  usual  form  including  an  electronic  part 

H  ,  a  pure  phonon  part  H  ,  and  an  interaction  terra  H  .  : 
e  ph  e*ph 
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H  =H  +H.+H 
l  e  ph  e • ph 


(2) 


It  is  useful,  following  Licciardello,  Stein  and  Haldane  (1980),  to  renor¬ 
malize  the  pure  terms  by  introducing  displaced  phonon  operator  d^  through 
the  relation 


di  =  bi  -  g  ni/w 

This  immediately  gives  the  form 


h.  =  a  ,  +  h  . 

i  el  ph 


where 


and 


H  ,  =  e.  n,  -  -r-  C  n.n.  +  Un. .n. . 
el  l  i  2  li  1+1+ 


«Ph  H  “  di  di 


tfy...-  _ ...  ... 
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(3) 


(4) 


(5) 


(6) 


and  where  we  have  defined  the  coupling  energy  C  =  2  g2/t»)  .  The 
stationary  solutions  are  described  by  quantum  numbers  n^,  which  give 
the  number  of  electrons  and  phonons,  respectively,  localized  at  the  site  i 
We  define  a  local  equilibrium  state  as  one  which  satisfies 


H  ,  n. ,v,>  =0 

ph1  i’  i  eq 


(7) 


or 


di  E  0 


It  is  the  purpose  of  this  note  to  examine  conditions  on  the  Hamiltonian 

parameters  under  which  the  state  In. ,v.>  is  metastable. 

'  i  i  eq 

We  consider  two  sites  i  and  j  whose  corresponding  states 

lni’Ui>eq  and  lnj’vj>eq  each  satisfy  equation  (7)  .  At  low 
temperatures,  we  may  roughly  define  stability  of  the  equilibrium  states  i 
and  j  with  respect  to  spontaneous  charge  transfer  between  them  by 


-3- 


requiring  the  energy  for  the  process  to  be  endothermic,  i.e. 


<n  v  H  n  v  >  +  <n  v  H.  n.v.> 

eq  i  i1  i1  i  i  eq  eq  j  j 1  j 1  3  3  eq 


<  cn^+An^^lH  |n  ♦An^^ 


+  <n .-An .v . I H . I n.-An . v .> 
3  3  3  3  3  3  3 


(8) 


where  An  =  ±1,±2  J  0  <  n .  .  +  An  <  2  .  It  is  instructive  to  consider 

-  1,3  ~ 

the  range  of  electronic  energies  e.  under  which  eq.  (8)  is  satisfied 
for  diamagnetic  n^  =  2  and  paramagnetic  states  n^  =  1  .  For  the  system 
in  equilibrium  defined  by  eq.  (7)  for  each  site,  the  energy  eigenvalue  is 
directly  given  through  equation  (5)  for  the  site  i  : 


E 


i 


en^  + 


U6 


2,n. 


(9) 


where  n^  is  the  number  of  electrons  at  the  site  i  and  6  is  the 

Kronicker  delta.  As  has  been  noted  (Anderson,  1975),  the  system  ground 

state  is  determined  by  the  occupation  set  {n^}  =  {0,2}  and  for 

convenience  we  define  the  Fermi  energy  e  by  the  relation 

r 

2ep  +  U  -  2C  =  0  (10) 

Thus,  in  the  ground  state,  sites  for  which  have  occupation  n^  =  0 

and  sites  for  which  e.  <  e_  are  occupied  by  n.  =  2  . 

1  r  I 

Infinite  Range  Model: 

We  define  a  metastable  state  i  as  an  equilibrium  state  i  (cf  eq.  7) 
for  which  with  n^  ^  0  or  with  n^  +  2  and  which  is 


i 


stable  according  to  eq.  (8)  for  all  j  .  The  latter  condition  defines  an 
infinite  range  model  which  we  consider  here.  This  assumption  has  been 
relaxed  in  other  works  (Licciardello,  Haldane,  Stein,  1980;  Licciardello 
1980b) . 


Metastable  Diamagnetism: 

We  consider  a  metastable  state  i  which  contains  n.  =  2  .  The 
bandwidth  for  stability  is  given  through  eq.  (8)  and  we  evaluate  the  case 
Uj  =  0  and  An^  =  -1  .  It  has  been  shown  that  the  two-electron 
emission  process  An^  =  +  2  is  always  more  severe  (i.e.  states  are  more 
stable  energetically  with  respect  to  this  process)  so  that  it  may  be 
ignored.  We  obtain,  using  eqs.  (8),  (4),  (5)  and  (6) 

-  Ej  <  2C  -  U  (11) 

as  the  condition  for  metastable  diamagnetism.  The  region  is  sketched  in 
Figure  1  where  the  stable  regions  ei  <  eF  and  ej  >  £f  *  and  the  Fermi 
point  ,  defined  through  eq.  (10),  are  also  indicated.  The  metastable 
states  occur  generally  in  the  shaded  triangle.  In  the  infinite  range  model 
which  assumes  non-zero  electronic  matrix  elements  between  all  i,j  ,  it  is 
easy  to  show  that  the  allowed  solution  corresponds  to  the  largest  rectangle, 
with  sides  parallel  to  the  axes,  which  may  be  inscribed  in  the  metastable 
region.  This  is  also  indicated  yielding  the  bandwidths 


for  n^  =  2  and 


C  -  U/2  <  ei  <  2C  -  U 
0  <  tj  <  C  -  U/2 


(12) 


(13) 


for  n^  ■  0 


in  agreement  with  the  results  of  Ref.  2.  It  is  interesting 
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Co  note  that  any  relaxation  of  the  range  of  the  interaction  increases 
the  range  of  metastability.  In  particular,  if  only  decays  to  the  ground 
state  are  permitted  (as  may  be  arranged  experimentally  for  a  system  with 
few  metastable  states)  the  bandwidths  are  increased  by  a  factor  of  2 
(indicated  by  the  lengths  of  the  sides  of  the  isosceles  triangle). 


Metastable  Paramagnetism: 

We  also  consider  metastable  states  occupied  by  a  single  electron 
which,  as  described  above  .have  no  stable  regions.  Thus  we  consider  a 
site  i  where  n^  =  1  and  consider  the  absorption  An  =  1  and  emission 
An  =  -1  of  a  single  particle  through  a  process  involving  the  site  j  . 
Again  using  the  criterion  (8)  to  define  the  region  of  metastability,  the 
range  of  energy  e ^  is  given  througi  the  relations 


nj  *  0 


ei  "  ej  -  C 


(14) 


n.  =  1 
3 


ei-£j*U 


(15) 


=2. 


VEiiC 


(16) 


where  we  have  employed  eqs.  (4),  (5),  and  (6).  The  paramagnetic  case  is 
sketched  in  Figure  2,  where  the  lines  defining  the  regions  satisfying  the 
inequalities  (14)  and  (16)  are  drawn.  The  vertical  boundaries  of  the 
shaded  region  where  paramagnetic  metastable  states  may  occur  are  defined 
through  the  limits  of  n^  =  2  states  (eq.(12)  ,  where  we  use  i  -+•  j) 
within  the  infinite  range  approximation  and  correspondingly  for  n^  =  0 
(eq.  (13)).  These  are  the  vertical  lines  at  =  2C  -  U  and  at 
Cj  *  0  ,  The  same  algorithm  applies  for  determining  the  allowed  paramagnetic 


sfates  which  is  the  largest  rectangle  with  sides  parallel  to  the  axes. 

This  region  is  indicated  as  (a)  for  the  infinite  range  model.  The  band¬ 
width  satisfies  6e^  <  2U  ,  as  is  seen  by  inspection,  first  predicted  by 
Licciardello  et  al.,  and  is  consistent  with  the  ineqality  (15)  without 
further  restriction. 

ignoring 

It  is  interesting  to  examine  the  consequences  of  Afast  particle 
exchange  between  paramagnetic  states,  i.e.  process  (15).  This  is  entirely 
realizable  at  densities  in  which  paramagnetism  has  been  induced  in  amorphous 
materials  up  to  now.  If  the  range  of  metastable  diamagnetic  states  is 
suppressed,  (indicated  in  rectangle  b  ,  in  Figure  2)  the  energy  range  of 
allowed  paramagnetic  states  increases  as  shown.  This  is  an  interesting 
prediction  and  may  be  tested  experimentally.  Ingeieral,  metastable  diamagnetic 
states  are  produced  thermally  (no  gap)  whereas  metastable  paramagnetic  states 
must  be  optically  induced.  Thus  materials  annealed  at  low  temperature  for 
long  times  should,  according  to  these  considerations,  be  susceptible  to  a 
larger  density  of  induced  paramagnetism. 
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Figure  1.  Phase  diagram  for  diamagnetic  metastable  states  with  energy 


and  occupation  n^  =  2  to  fast  emission  of  an  electron 
to  a  site  g  with  n^  =  0  and  energy  .  The  rectangle 
defines  the  metastable  region  when  the  process  i  -*■  j  is 
permitted  for  all  i,j  . 

Figure  2.  Phase  diagram  for  paramagnetic  metastable  states  with  ni  =  1 
and  energy  .  The  shaded  region  defines  the  region  of 
stability  for  such  states  to  fast  emission  of  absorption  of  a 
single  electron  from  metastable  states  with  energy  . 

The  rectangular  area  (a)  defines  the  metastable  region  when 
the  process  i  -*■  j  is  allowed  for  all  i,j  .  The  region  (b) 
ignores  processes  involving  only  paramagnetic  states  and 
considers  processes  involving  a  narrower  region  of  metastable 

states  with  n.  =  0,2  . 
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